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The continual model of the nonuniform magnetism in thin films and wires of a diluted 
magnetic semiconductor is considered with taking into account the finite spin polarization of 
carriers responsible for the indirect interaction of magnetic impurities (e.g. via RKKY mechanism). 
Spatial distributions (across the film thickness or the wire radius) of the magnetizaton and carrier 
concentrations of different spin orientations, as well as the temperature dependence of the average 
magnetization are determined as the solution of the nonlinear integral equation. 

The indirect interaction of magnetic impurities, e.g. of the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) type, is believed to be one of the basic mechanisms of magnetic ordering in systems 
with free carriers of a high concentration (metals and degenerate semiconductors). Since the 
most of potentially interesting electron devices are of nanometer sizes, it is of interest to see 
how magnetic properties of the relevant systems depend on their finite sizes [HE]- The typical 
example is the thin film of the diluted magnetic semiconductor (e.g., Gai-^Mn^As) with the 
thickness on the order of some tens of the lattice periods which are considered in the present 
paper. In such systems, the mobility of charge carriers (holes) is low (on the order of 10 
cm 2 /V-s [3]) and, hence, their mean free path is also short. Then the collisional broadening 
of the hole energy levels is so high that the system appears to be effectively three dimensional 
one, and there is no need to take into account effects of the size quantization [4]. 

The magnetization of a thin enough film of the diluted magnetic semiconductor (whose 
thickness L is comparable with the characteristic length £ of the magnetic impurities' indirect 
interaction, see below) could be essentially nonuniform, at least, in the direction perpendicular 
to the film plane. It is convenient to characterize such a nonuniform magnetization, depending 
on coordinates, by the reduced local magnetization j(r) = M{v)/M s (M s is the saturation 
magnetization). Non-zero local magnetization (0 < j < 1 ) of Mn atoms with the spin 
SWi = 5/2 leads to the nonuniform local spin polarization of holes reflecting in the fact 
that the local concentration p~(r) of holes with spins, antiparallel to the local magnetization, 
exceeds the concentration p + (r) of holes with the opposite spin direction. At that, the local 
hole polarization £(r) = [p~(r) — p + (r)]/p(r) does not equal zero (p(r) = p~(r) +p + (r) is the 
total hole concentration). Formally, the hole polarization could be related with the effective 
spin-depending magnetic potential [5] 

Vmag(r) = xN a 3 J pd aSunj(r), (1) 

influencing on the holes. Here Ao ~ 8 • 10 21 cm -3 is the concentration of Ga sites in the GaAs 
lattice, x is the fraction of such sites populated by Mn atoms, a being the lattice period, and 
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a = ±1/2 is the hole spin, J p d = 1.2 eV is the exchange energy between mobile holes and 
localized (^-electrons of Mn atoms \6\. According to (DQ), major holes (with the preferred spin 
orientation) are accumulated in regions of the high magnetization, and, contrary, minor ones 
are pushed out in regions with the low magnetization. 

Such a spin separation of carriers in the space could be, in principle, accompanied by 
modifying the spatial distribution of holes resulting in the violation of the local charge neutrality. 
To see if that possible in the considered case, notice that the characteristic displacement of the 
hole charge "centroid" over a distance 5L leads to the electric potential equal to 
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V* ~ ^-{6L)\ (2) 

K 

where k ~ 10 is the semiconductor dielectric constant. From (DO), ([2j) one finds Ki/ 1 Knag I ~ 
e 2 L 2 /J pd . Herefrom, it follows V^i/|V^n a g| 1 at 5L > 10 A. This suggest that for films of 
the thickness L ^> 10 A one could neglect the violation of the electric neutrality and hold 
p(r) = p~(r) + p + (r) = Const. In the present paper, we shall not go beyond considering such a 
case. 

Let the Fermi energy of charge carriers in the absence of the magnetization (j = 0) be e F . 
Magnetic potential (pQ), arising with the origination of the magnetization, leads to splitting the 
hole band into two spin sub-bands with effective Fermi energies (measured from bottoms of 
those sub-bands) e F + Kn ag (r) and e F — V mag (r) which correspond to following Fermi momenta 

k$(r) = kp(l±^^) 1/ \ (3) 



e F 



where kp = {2m*e F /h ) ' , m* w 0.5m is the effective hole mass (at V mSLg (r)/ep > 1, 
k F (r) = 0). 

Carrier concentrations in the sub-bands are determined by the set of simple equations 

p-(r)(x[e F + V mag (r)f 2 , p+(r) oc [e F - \/ mag (r) f 2 , p~(r) +p + (r) = p, (4) 
from which it follows 

P-W = ?[! + £(')], P + W = ?[l-eW], (5) 



where 



2 " v /J ' ' v ' 2 



, , u(r) — 1 



is the hole spin polarization, u(r) = {[1 + V mag (r)/e F ]/[l - V mag (r)/e F }} 3/2 (if V mag (r)/e F > 1, 
then p~ = p, p + = 0, £ = 1). 

With the uniform magnetization (j = Const), momenta kp do not depend on coordinates 
and the expression for the energy w(p) of the indirect interaction of two magnetic atoms, 
separated by a distance p from each other, taking into account the spin splitting of the hole 
band could be written in the form 

w(p) = -~I $(p)exp(-p/e), (7) 
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where, as an example, the RKKY range function reads [7] 

$(p) = Q [0+(p)cos0 + (p) -sin0+(p) + 0-(p)cos0 _ (p) -sinr(p)], (8) 

^^f^*)' 9± ^ = 2k rP- ( 9 ) 

Exponential factor in (JTJ) allows for the finite length £ of the hole spin relaxation [8] (in the 
simplest case, coinciding with their mean free path). 

With increasing magnetic potential (JTJ) , the intensity of the interaction (jT]) is changed. One 
could judge that from the behaviour of the function <3>(p) at distances close to the mean distance 
p between impurities. Spatial dependencies of the RKKY-function $(p) for various values of 
the ratio V mag /eF are shown in Fig. 1. It is clear that in the inter-impurity distance range of 
p = 1.5a — 2a, corresponding to the impurity concentration x = 0.05 — 0.1, the interaction 
energy (at a fixed Ef value) decreases noticeably with increasing Vm a g. What that means is the 
spin splitting of the hole band has to result, in the end, in decreasing the system magnetization. 

The generalization of the function (fSj) over the case of the non-uniform magnetization could 
be done with replacing phases ± (p) by their average values 

p p 
0* (p) = 2 J k$(s)ds = 2k F J ^l±Aj(s)ds, (10) 



where A = xN a 3 Jpd&Syin/tF, and the integration is performed along the line connecting the 
impurities. Then the function $(/?) becomes the functional of the spatially non-uniform magne- 
tization: 



, 4 

a 



P 



${p)exp{-p/e)^F[Aj(T),p] = 
[6* + {p) cos0;(p) - sin0;(p) + 6*_{p) cos6*_{p) - sm6*_{p)]e~ p/e . (11) 



Under low magnetization (J <C 1), the magnetic potential (DD) is also low, and the expression 
( fTTi ) goes to the standard form. 

The finite carrier polarization is important in the only case when the magnetic potential 
and Fermi energy are comparable, that is at Aj ~ 1. Let us estimate the parameter A for 
the magnetic semiconductor Gai_ x Mn x As. Due to the compensation, the hole concentration 
p is always lower than the concentration xN of Mn atoms (acceptors) positioned in Ga-sites. 
Nevertheless, p/(xN ) > 0.3 at x = 0.05 [9]. The estimate for that case leads to A > 1. If, 
in addition, the local magnetization in some regions of the system is close to the saturation 
(j ~ 1), then Aj ~ 1 and taking account of the finite carrier polarization could be important. 

To test the validity of that conclusion, we have carried out calculating the spatial distribution 
of the magnetization in a thin magnetic semiconductor film where the magnetization is parallel 
to the film surface. In doing so, we have employed the approach elaborated early [4]. 

The energy H^rkky of the indirect interaction of a given spin Si with its surraundings is 
defined by the relation Wrkky = ^2jLi w (Pij)' The distance p^ could not be less than the 
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distance r min between the neighbor lattice sites, accessible for magnetic impurities (for the 
diluted semiconductor Gai-zMn^As, r min = a/ a/2 where a « 5A being the lattice constant). 



In the continual approximation, the sum could be replaced by the integral 
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W RKKY (r) = -^/o / F[Aj{v'), |r - r'\}j(r')d 3 r' 



where n M is the concentration of magnetic Mn-ions, and the integration is expanded over the 
volume occupied by impurities. 

Contrary to the infinite system, the value of that integral depends on the system geometry 
and the position of a chosen point. Below, we consider the magnetic interaction in the film 
of the thickness L (in the z-axis direction) which is so thin that the local magnetization j(h) 
within the film is always parallel to the surface and depends on the given point distance h from 
the median film plane onl)lE In that case, it follows from (112]) 



L/2 



RKKY 



(h) 



\ 



F[Aj(z),p}pdp 



Z=-L/2 \p m in(z,h) 



j{z)dz, 



(13) 



/ 



where 



Pmin(z, h) = Max[ I h - z\, rmin]. (14) 

The relation (Tl3l . defining the local energy WrkkyCO in all film points, is non-local one: 
that energy is the functional of the local magnetization j(p) and specified by all film points. 

The self-consisted equation j(h) = ta.iah[W-[iKKY{h)/kT], determining the local magnetization, 
could be now written as follows 

L/2 ( oo \ 



j(h) = tanh 



kT 



-L/2 



F[Aj(z),p}pdp 



j(z)dz 



J 



(15) 



where phases appearing in fTTT]) are defined by relations 



9* ± (p, z, h) = 2k F J y/l±Aj[h-(f//p)(h-z)] dp'. 
Introducing the designation 



(16) 



- nn^ J F[Aj(z),p]pdp = K(z, h) 

Pmin(z, h) 

for the "internal" integral in Eq. (fT5] ). one could write it in the form 

L/2 
1 f 

j(h) = tanh — / K(z,h)j(z)dz 

z=-L/2 



(17) 



(18) 



2 That means we consider systems with weak enough surface magnetic anisotropy [10] 
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where r = kT / Iq is the reduced temperature. 

The non-linear integral equation (fT8l) determines the spatial distribution of the magnetization 
in the considered system at a given temperature. It differs from the equation considered early|4j 
in that its kernel is the functional of the non-uniform magnetization (due to the finite carrier 
spin polarization). The solution of the obtained equation has been found, as in [4j, by the 
successive approximations method. 

Spatial distributions of the magnetization in the film of the thickness L = 20a at various 
temperatures r are shown in Fig. 2. The chosen value k F a = 1 corresponds to the hole 
concentration p rs 2 • 10 20 cm -3 , characteristic to Gai_2;Mn x .As with x = 0.05. It is clear 
that taking account of the band splitting, leading to the finite carrier polarization, decreases 
noticeably, as might be expected, the system magnetization. Near the film surface the magne- 
tization drops quickly, and the thickness of the relevant subsurface regions is almost linearly 
increase with the temperature (cf. Fig. 2). 

The magnetization of "thick" films (L — ► oo) is nearly uniform (except for "thin" near- 
surface regions where it is twice as low comparing to the bulk) . The effect of the band splitting 
evaluated for that case by the ratio j(A — 1,L — > oo)/j(A = 0,L — > oo), is illustrated by 
Fig. 3, where its temperature dependence is shown. It could be seen that at r > 0.2rc the effect 
is very notable. 

It is convenient to characterize the non-uniformly magnetized film by its average magnetization 

L/2 

U) = j; I j(z)dz- (19) 

-L/2 

Temperature dependencies of that quantity for the film of the thickness L = 20a are represented 
in Fig. 4. They allow to find the Curie temperature Tq, defined as the temperature at which 
0) ~~ * 0- Near the Curie temperature, j < 1 and, hence, taking into account the finite spin 
polarization does not change Tc value. That is immediately seen from Fig. 4, where both 
dependencies (j)(r) (for A = and A = 1) give the same value tc ~ 1. 

In the framework of the considered scheme, it is easy to find spatial distributions of the 
carrier spin polarization, as well as those of major and minor carrier concentrations. Correspon- 
ding dependencies, obtained with the help of Eqs. flE},®, are shown in Fig. 4. It follows, 
therefrom, that the hole polarization £ in the central part of the film is more higher than at 
the periphery. That is connected with the spatial separation of carriers with different spin 
polarization: major carriers concentrate away from the film surface, and minor ones - in 
peripheral regions. (In those regions, the violation of the electric neutrality is possible, so the 
relevant results should be considered as qualitative ones.) 

At T — > Tq Eq. (fl~8l) is simplified, turning into the homogeneous linear integral equation: 

L/2 

j{h) = - [ K(z,h)j(z)dz, (20) 

T"C J 

z=-L/2 

from which one could see that the Curie temperature tq is nothing than the eigenvalue of that 
equation kernel. Due to the symmetry of the kernel (K(z,h) = K(h,z), see. (fl~7l) ) such an 
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eigenvalue is always exists [TT]. To determine it, any one of a number of known numerical 
methods [12 ] is appropriate. However, in practice it is more simply to find an approximate tq 
value by the above described approach. 

In the same way, magnetic properties of a thin wire could be considered. The generalization 
of the relevant integral equation [IB], for the case of the local magnetization being everywhere 
oriented along the wire axis and depending on distance h to that axis only, reads 



j(h) = tanh 

where R is the wire radius, 



T 

L o 



R 

K(r, h) j(r)dr 



(21) 




K(r, h) = -n,r \ \ / ===== 7^ I <V, ( 22 ) 

'/r — r A — h z + 2rn cos cp 

Pmin = Pmm(r, h, 4>) = Max[(r 2 + h 2 — 2r/icos0) 1 / 2 ,r min ]. In that case, phases, appearing in 
Eq. ( TTT1) . are defined by the relations 



9± = 2k } 



\j 1 ± Aj[{[h(l - p'/p) - r(p'/p)] 2 + 2r%'/p)(l - p'/p)(l - cos0)} 1/2 ] dp'. (23) 



In Fig. 5, spatial distributions of the local magnetization are displayed for the wire of the 
diameter 2R = 20a where at r ~ 1 the paraxial region occurs to be magnetized only. As 
well as for a film, the splitting of the carrier band results in the essential lowering the system 
magnetization. E.g., at r ~ 0.1 the magnetization in points, spaced at the distance h = R/2 
from the axis, decreases by 2.5 times because of that effect. 

The strong non-uniformity of the local magnetization results in that temperature dependen- 
cies of the average wire magnetization, shown in Fig. 6, have unusual concave shape. In the 
inset of Fig. 6, spatial dependencies of the concentrations p~, p + of major and minor carriers 
along with their polarization are displayed. One could see that the noticeable hole polarization 
occurs in the paraxial wire region only and (at r = 0.1) reaches ~ 80%. 

In conclusion, we have considered the continual model of the non-uniform magnetism in 
thin films and wires of diluted magnetic semiconductor under a finite spin polarization of 
carriers responsible for the indirect (RKKY) interaction of magnetic impurities. We succeeded 
in obtaining spatial distributions of the magnetization described by the non-linear integral 
equation with the help of quickly converging iterative procedure. They occurs to be essentially 
non-uniform and strongly dependent on splitting the band of carriers responsible for the 
indirect interaction between magnetic impurities. Results could be used to describe properties 
of nanosized systems of diluted magnetic semiconductors. 
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Captions 



Fig. 1. Range functions <&(p) of the inter-impurity interaction energy at various values 
of the ratio V majg /eF- Accepted: k F a = 1. 

Fig. 2. Spatial distributions of the local magnetization in the film of the thickness L = 20a 
at various temperatures r with [A — 1, right panel) and without (A = 0, left panel) taking 
account of the band splitting. Accepted: £ = 3a, kpa = 1, 47m M = 1. Areas of the reduced 
peripheric magnetization are marked by dashed lines. 

Fig. 3. Temperature dependence of the band splitting effect leading to the noticable reducing 
the magnetization of the "thick" (L — > oo) film. Accepted: £ = 3a, kpa = 1, Airn^ = 1. 

Fig. 4. Temperature dependencies of the average magnetization (j) in the film of the 
thickness L = 20a with (A = 1, lower curve) and without (A = 0, upper curve) taking account 
of the band splitting. Insert: spatial distributions of the spin polarization £, and concentrations 
of major (p~) and minor (p + ) carriers at A = 1, r = 0.8. Accepted: I = 3a, kpa = 1, 47m M = 1. 

Fig. 5. Spatial distributions of the local magnetization in the wire of the diameter 2R = 20a 
at temperatures r = 1 and r = 0.1 with (A = 1, solid curves) and without (A = 0, dashed 
curves) taking account of the finite carrier spin polarization. Accepted: £ = 3a, kpa = 1, 
47™^ = 1. 

Fig. 6. Temperature dependencies of the average magnetization (j) in the wire of the 
diameter 2R = 20a with (A = 1, lower curve) and without (A = 0, upper curve) taking account 
of the band splitting. Insert: spatial distributions of the spin polarization £, and concentrations 
of major (p~) and minor (p + ) carriers at A = 1, r = 1. Accepted: £ = 3a, k F a = 1, 47m M = 1. 
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